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Abstract 

We present a detailed analysis of XMM archival data of two Ultraluminous X-ray 
Sources (ULXs) in the nearby spiral galaxy NGC 1313: NGC 1313 X-l and X-2. The 
spectral continuum of these sources was modeled with a soft thermal component plus 
a power-law. If the soft component originates from an accretion disk, the inferred 
mass of the compact remnant is }t IOOMq, making it an Intermediate Mass Black 
Hole (IMBH). A detailed analysis of the residuals of the XMM EPIC-pn spectrum 
shows some evidence for the presence of an Oxygen emission line in NGC 1313 X-l. 
The simultaneous presence of an excess in emission, although at a much reduced 
significance level, at different energies in the X-ray spectra of NGC 1313 X-l and 
X-2 is suggestive of typical emission lines from young supernova remnants. An 
optical counterpart for NGC 1313 X-2 was also identified. On an ESO 3.6 m image, 
the Chandra error box embraces a R ~ 23 mag stellar-like object and excludes a 
previously proposed optical counterpart. 

Key words: Black holes, X-rays: individual (NGC 1313 X-l, NGC 1313 X-2/MS 
0317.7-6647), X-rays: binaries, X-rays: galaxies 



1 Introduction 

Point-like, off-nuclear X-ray sources with luminosities in excess of the Edding- 
ton limit for one solar mass are increasingly discovered in nearby galaxies. 
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Despite the growing body of observational data, the mechanism which pow- 
ers ULXs is still under debate. In several variability on timescales 
of months/years has been detected, hinting towards the presence of a com- 
pact object. Were these sources X-ray binaries in the host galaxy, assuming 
Eddington-limited accretion, masses in the range Mbh ~ 100 — 1OOOM are 
inferred for the compact object from the observed flux. IMBH may have orig- 
inated from the collapse of a massive star formed in a low metallicity environ- 
ment, or through mergi ng of massive stars (or lower mass BHs) in a cluster 
fjMiller fc Colbe rt. 2004). The BH mass estimates derived from considerations 
involving the Eddington limit are, however, questionable because these sources 
need not to be spherically symmetric, nor stationary. It has been proposed that 
many of the ULX proper t ies ca n be explained assuming that they do not emit 
isotropi cally (iKing et al!l200lh or are dominated by emission from a relativis- 
ts jet ( Kaaret et al. . 2003j ). In this case, they may harbor stellar mass BHs 
and may be similar to Galactic microquasars. 

Particularly interesting is the case of the three ULXs hosted in the nearby 
SBc galaxy NGC 1313. One of them is known to b e associated with the inter- 
acting supernova SN 1978K ( Schleeel et al. . 2000h . The other two have been 
extensively studied with several X-ray telescopes and, to date, provide some of 
the best evidence for the presence of a soft component in the X-ray spectrum 



of UL Xs. Moreover, as discussed in iMiller et al.l (J2004J) and IZampieri et al. 
(2004 ), hereafter Z04, their continuum subtracted XMM EPIC-pn spectra 
show significant residuals especially in the 0.5-3 keV range. 



2 X-ray Data 



XMM-Newton observed NGC 1313 on October 17, 2000 for a total of ~ 42 
Ks. The field was centered on the galaxy nucleus and contains all the three 
ULXs. The three XMM EPIC cameras operated in Prime Full Window mode 
with the medium filter. The analysis reported here follows that presented in 
Z04. We consider only the EPIC-pn spectrum because it has twice as many 
counts as each single MOS instrument. The EPIC-pn spectrum was directly 
extracted from the observation data file because the automatic pipeline pro- 
cessing failed to produce an event list. Data screening, region selection and 
event extraction were performed with the standard software XMM- S AS v 
6.0.0. In order to eliminate the possible contamination of solar flares (present 
during the observation), event files were filtered using the good time intervals 
when the total off-source count rate above 10 keV is less than 1 counts s _1 . 
This leave ~ 21 ks of useful time with an average count rate of 0.73 counts 
s- 1 for NGC 1313 X-l and 0.24 counts s" 1 for NGC 1313 X-2. 

The best fit of the XMM EPIC-pn continuum for both sources is obtained with 
an absorbed soft, thermal component (a multicolor disk blackbody, MCD) 
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plus a power-law. In comparison with previous RO SAT and ASCA data the 
statistics is sufficiently good that two components models provide a significant 
improvement over single component ones (see Z04). The statistical improve- 
ment obtained adding a MCD component to a power-law model is significant 
above the 4.5a level. The resulting best fitting parameters for NGC 1313 X-l 
and X-2 are reported in Table 1. The inferred inner disk temperature is much 
lower th an that obta i ned fo r single component thermal models, as already 
found bv iMiller et al. The X-ray luminosity in the 0.2-10 keV range, 



assuming isotropic emission and a distance of the host galaxy ~ 3.7 Mpc, is 
~ (1.4 ± 0.2) x 10 40 erg s" 1 for NGC 1313 X-l and ~ (6.0 ± 0.5) x 10 39 erg 
s _1 for NGC 1313 X-2. If at maximum the source radiates at the Eddington 
limit, the BH mass is ~ 120 M for X-l and ~ 50 M for X-2. Sub-Eddington 
accretion would imply an even larger mass. 

Residuals in the EPIC-pn spectra of NGC 1313 X-l and X-2 suggest the 
presence of some emission lines in both sources. In order to test for the presence 
of spectral features, we added, one by one, gaussian components at different 
energies where residuals show evidence of some excess in emission. 

For NGC 1313 X-l, the most significant residual is at around 0.6 keV, identified 
with a high ionization (He-like) Oxygen line (see Figure 1). Adding a gaussian 
component, the improvement of the fit with respect to an absorbed MCD+PL 
model is significant at the 3.5 a level. The line energy is 0.58 keV, while the 
line width was frozen in the fit (0.01 keV). Residuals show also an excess of 
emission at energies of 1.8 and 4.7 keV (see Figure 1). The first is the typical 
energy of highly ionized Silicon, while the second is not readily identified. 
Adding other two gaussian components for fitting these residuals does not 
give any further statistical improvement in the fit. 

A similar analysis was performed also on the EPIC-pn spectrum of NGC 
1313 X-2, following the same approach outlined above. Also in this case the 
spectrum shows some residuals at an energy of 0.6 keV. In order to test the in- 
flu ence of ISM absorp tion models (due to the neutral Oxygen absorption edge, 



sec 



Miller et al.ll2004h on the detection of an O line in NGC 1313 X-l, we re- 
peated our analys is using the Tubingen-Boulder absorption model (TBABS; 
Wilms et al. 2000L Table 1), that includes a treatment of gas-phase ISM with 



revised photoionization cross sections and revised abundances (plus the contri- 
bution of grain-phase ISM and molecules). The improvement obtained adding 
a gaussian component to a TBABS+MCD+PL continuum is at the ~ 3.6 a 
level. The line energy is 0.6 keV, while the line width was frozen in the fit (0.01 
keV). The significance of the line remains fairly high even performing a fit in 
a restricted energy region (0.3-1.0 keV) around the centroid of the gaussian. 
Fitting the continuum with an absorbed MCD+PL model, the addition of a 
gaussian component with fixed width (0.01 keV) gives an improvement at the 
~ 3 a level. 
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Additional information on the line properties were obtained analyzing the line 
profiles. After freezing the best fitting continuum of Table 1 and removing the 
gaussian components, residuals (normalized to the continuum) were fitted with 
a constant (equal to unity) plus a gaussian profile, obtaining the lines centroid 
energies and equivalent widths. The energies are consistent (within the errors) 
with those derived from the spectral fits, while the equivalent widths are all 
below 100 eV. A thorough analysis of these data and new XMM observations 
of NGC 1313 X-l and X-2 will be presented in Mucciarelli et al. (2005, in 
preparation) . 



3 The optical counterpart of NGC 1313 X-2 



Optical images of the field of NGC 1313 X-2 in the R-band were taken on 
16 January 2002 with the 3.6 m telescope of the European Southern Ob- 
servatory (ESO) at La Silla (Chile). Four images were obtained for a total 
exposure time of 1320 s (see Z04 for details). The accuracy of the astrometric 
calibration, performed with GSC2 ESO field stars, is 0.3" (1-er). The X-ray 
position of NCG 1313 X-2 was obtained from a 2002 Chandra pointing with an 
accuracy of 0.7" (Z04). The astrometric calibration of the optical and Chan- 
dra images was che cked using the very accurate radio position of SN 1978K 
fjRvder et al.lll993l Z04). Our Chandra posi tion of NGC 1313 X -2 is shown 
in Fig ure 2, together with the ROSAT HRI (|Schlegel et all l2000h and XMM 
EPIC (|Miller et all [2003) error boxes, overlaid on our ESO image. All mea- 
surements are consistent within l-a. Object C (R ~ 23) is inside the Chandra 
error box and its position coincides within l-a with that of NGC 1313 X-2, 
making it a likely counterpart and ruling out the previous proposed counter- 
part, object A. From the maximum absorbed X-ray flux of NGC 1313 X-2 
(fx ~ 2 x 10~ 12 erg cm -2 s _1 ) and optical magnitude of object C, we estimate 
fx/ fo P t £ 1000. 



4 Discussion 



The X-ray spectral parameters, in particular the temperature of the MCD fit 
(Tmcd), can be used to estimate the BH mass. From the effective temper- 
ature of a standard accretion disk, under the assumption that T MC d repre- 
sents an estimate of the maximum disk temperature, Z04 infer Mbh/Mq = 
(Mc 2 1 ' LEdd)f A (oiTMCD/^-^ x 10 7 K) -4 (where / ~ 1.5 is a color correction 
factor). The temperature obtained from the two-components fit to the XMM 
spectrum (Table 1) implies M BH rj (130/90)/ 4 M for NGC 1313 X-l/X-2 
(assuming a ~ 2 and Eddington limited accretion). 
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Up to now evidence of em ission lines in UL Xs has been reported in M82 X-l 
" 2003lh NGC 4559 X-10 ((Cropper et all l2004h and M101 



f Strohmaver et al.. 


2003IL 


ULX-1 (Koner et al. 


. 200f# 



spectrum of NGC 1313 X-l appears to be significant at above 3.5 a while the 
statistical evidence of other features is at most marginal (< 2 a). The simplest 
interpretation is that we are observing typical emission features of intermediate 
mass elements left over after the explosion of one or more supernovae belonging 
to the stellar association of massive stars where the two ULXs are probably 
embedded. 



Concerning NGC 1313 X-2, the luminosity inferred from the apparent mag- 
nitude of object C is consistent with a ~ 2OM main sequence star or a 
~ 15 - 2OM evolved OB supergiant making NGC 1313 X-2 a High-Mass 
X-ray Binary. This picture seems to be confirmed also by the study of the 
environment of NGC 1313 X-2 (see Z04) and by the association o f both ULXs 
with extended optical emission nebulae (|Pakull fc Mirioni . 20021 ). 
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Table 1 

Parameters of the fit of the XMM-EPIC observation of NGC 1313 X-l and X-2 



Model 


1() 21 ,.)« - 


r 


kT [keV] 


Parameters 


X? cd (dof) 


NGC 1313 X-l 


WABS(MCD+PL) 

TBABS(MCD+PL) 

TBABS(MCD+PL+GAUSS) 


n 4 o + 0.60 
•*- 48 -0.25 

4 21+ 64 
-0.32 

4 48+ - 73 
4 -^°-0.31 


i 7C+O.O8 
1 69+ 07 

i 79 +0.06 


19+ 001 
19+ 001 
17+ - 01 


E = 0.60to.o2 ke V 


1.12(301) 
1.10(301) 
1.04(299) 


NGC 1313 X-2 


WABS(MCD+PL) 


3 26+ 92 
°- zo -0.25 


n 17+O.I5 
z - x '-0.09 


20+ ' 05 




1.21(109) 




Enerqy [kevl 

Fig. 1. XMMEPIC-pn spectrum of NGC 1313 X-l: absorbed MCD (dotted line)+PL 
(dashed line) model with three gaussian components (thin solid lines) at energies 
0.59, 1.8 and 4.7 keV. 
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Fig. 2. ESO 3.6m i2-band image of the field of NGC 1313 X-2. The circles show 
the ROSAT HRI, XMM EPIC and Chandra ACIS-S positions. The estimated 90% 
confidence radii are 6", 4" and 1.4" respectively. Labels A and C mark the old and 
new proposed optical counterparts. The insert at the bottomdeft shows a larger 
portion of the image with the position of the X-ray source (cross). 
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